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^ i' Abstract 

^^ - A collection is made of presently unexplained phenomena within our Solar system and in the 

universe. These phenomena are (i) the Pioneer anomaly, (ii) the flyby anomaly, (iii) the increase 
of the Astronomical Unit, (iv) the quadrupole and octupole anomaly, and (v) Dark Energy and 
(vi) Dark Matter. A new data analysis of the complete set of Pioneer data is announced in 
order to search for systematic effects or to confirm the unexplained acceleration. We also review 
r»o ' the mysterious flyby anomaly where the velocities of spacecraft after Earth swing-bys are larger 

l/^ , than expected. We emphasize the scientific aspects of this anomaly and propose systematic and 

Cp ■ continuous observations and studies at the occasion of future flybys. Further anomalies within 

the Solar system are the increase of the Astronomical Unit and the quadrupole and octupole 
anomaly. We briefly mention Dark Matter and Dark Energy since in some cases a relation 
f~j) ', between them and the Solar system anomalies have been speculated. 

I . 1 Introduction 

bX). 

Progress in physics always has been stimulated by observations which could not been explained 
within the presently standard physical theories. In the late 19th century the observations and 
^ ' experiments by Bradley as well as Airy who both observed aberration of distant starlight, of Fizeau 
j^ ■ who observed a dragging of light in moving media not compatible with theory at that time, and 
finally the experiment of Michelson and Morley who showed that the failure of the application of the 
non-relativistic mechanically laws to light propagation. All these effects which could not be made 
compatible with non-relativistic physics without introducing various unnatural elements into the 
theory, culminated into the invention of Special Relativity. Then, the theoretical incompatibility 
of Newtonian gravity with Special Relativity as well as the since long observed perihelion shift 
of Mercury which first has been attributed to systematic errors or solely to be due to the Solar 
quadrupole moment, lead to the formulation of General Relativity. Later on, the experimental 
study of atomic spectra which could not be explained using the laws of classical mechanics first led 
to Bohr's atomic model and, subsequently, to the the various formulations of quantum mechanics. 
The situation of gravitational physics today bears many similarities. At first, the theoretical 
inconsistency of quantum mechanics and General Relativity makes a new theory combining these 
two universal theories necessary. Furthermore, there are observations which at least until now and 
after many years of studies, have not yet found any convincing explanation. These observations are 
(i) dark energy which is necessary - under the assumption of the validity of Einstein's equations 
- to describe the accelerated expansion of the universe and (ii) dark matter which - again under 



the assumption of General Relativity - is necessary to account for the galactic rotation curves, for 
observed gravitational lensing of light, and for the structure formation in the early universe. Of a 
slightly weaker observational basis is (iii) the Pioneer anomaly, an unexplained constant accelera- 
tion of the Pioneer 10 and 11 spacecrafts, (iv) the flyby anomaly, an unexplained increase of the 
velocity of a series of spacecrafts after Earth Gravity Assists, (v) the recently realized increase of 
the Astronomical Unit defined by the distance of the planets from the Sun by approximately 10 m 
per century, and (vi) the quadrupole and octupole anomaly which describes the correlation of the 
low / contributions of the Cosmic Microwave Background to the orientation of the Solar system. 

These six phenomena, including Dark Energy and Dark Matter which at this stage are nothing 
more than a synonym for these observations, had neither found any convincing interpretation or 
solution nor culminated into a finally convincing theory. Lacking any explanation until now, these 
phenomena have the potential to be of importance for a new physics. 

In this paper we describe all these unexplained observations, state the open questions, and 
suggest new observations and new missions in order to obtain better data for a better analysis of 
these phenomena. 

2 Dark matter 

Dark matter has been introduced in order to " explain" the gravitational field needed for the galactic 
rotation curves, the gravitational lensing of galaxies, and the formation of structures in our universe 
[1] . It also appears in the spectral decomposition of the cosmic microwave background radiation [2] . 
Dark matter is needed if one assumes Einsteins field equations to be valid. However, there is no 
single observational hint at particles which could make up this dark matter. As a consequence, there 
are attempts to describe the same effects by a modification [3] of the gravitational field equations, 
e.g. of Yukawa form, or by a modification of the dynamics of particles, like the MOND ansatz [4, 5]^, 
recently formulated in a relativistic frame [7]. Due to the lack of direct detection of Dark Matter 
particles, all those attempts are on the same footing. 

3 Dark energy 

Similarly, recent observations of type la supernovae indicate that the expansion of the universe is 
accelerating and that 75% of the total energy density consist of a dark energy component with 
negative pressure [8, 9]. Furthermore WMAP measurements of the cosmic microwave background 
[10, 11], the galaxy power spectrum and the Lyman-alpha forest data lines [12, 13, 14] also indicate 
- when compared with standard cosmological models - the existence of the mysterious Dark Energy 
that leads to the acceleration of the universe, rather than a modification of the basic laws of gravita- 
tion [15]. However, also in this case there are attempts to give an explanation in terms of modified 
field equations, see, e.g., [16]. Recently it has been claimed that dark energy or, equivalently, the 
observed acceleration of the universe can be explained by inhomogeneous cosmological models, such 
as the spherically-symmetric Lemaitre-Tolman-Bondi model, see, e.g., [17, 18, 19]. 

^In a nice short paper M. Veltman [6] speculates how astronomers may build up laws of gravity by observing gravity 
on larger scales (scale of galaxies and of the universe) and compares that with a quantum field theory approach. 



4 The Pioneer anomaly 

The Pioneer anomaly is an anomalous unexplained acceleration of the Pioneer 10 and 11 spacecraft 
of 

apioncer = (8.74 ± 1.33) • 10"^° m/s2 (1) 

toward the Sun [20, 21]. This acceleration seems to have been turned on after the last flyby at 
Saturn and stayed constant within an 3% range. 

4.1 The observation 

The principle of observation was 2-way Doppler tracking: A sender on the Earth emits a signal of 
frequency fo which is "seen" by the spacecraft as frequency 

v' = , ^ f 1 - - Vo • (2) 

The spacecraft sends this frequency back (with a slight offset what, however, will not affect the 
principle of measurement), so that the receiver on Earth observes the frequency 

-" = ^^=(1--)-'. (3) 

Vl - v^/c^ V c) ^ ^ 

The comparison of the sent and received frequency gives the velocity of the spacecraft 

. = ^^^^^ = -2:^--2^. (4) 

Uq 1 + V/C C 

This measured frequency can be compared with the frequency obtained from the calculated orbit 
given by the gravitational field inside the Solar system together with all kinds of modeling needed 
(see below). 

The outcome the observation was a continuous drift between the observed and calculated fre- 
quency shift 

^^^^^^^^^^^^ = (3.84 ± 0.01) • 10-18 s-^ (5) 

This corresponds to a continuous drift in the velocity of the spacecraft or, equivalently, in a constant 
acceleration (1). 

4.2 Orbit determination 

Since the observations are made with tracking stations on the moving Earth observing the frequency 
of signals, the orbit determination consists of five segments: 

• Model of gravitational forces 

• Model of external non-gravitational forces 

• Model of internal (spacecraft) non- 
gravitational forces 

• Model of observation stations 

• Model of signal propagation 

• Codes 

We just mention the main aspects of this scheme. Most of this can be found in [21]. 



Orbit and velocity determination 



Gravitational forces The calculation of the orbits have been performed using relativistic equa- 
tions of motion for celestial bodies including order v^/c^ 

• The relativistic gravitational accelerations (EIH model) include the Sun, the Moon, and 
the 9 planets as point masses in an isotropic PPN A^-body metric. 

• Newtonian gravity from large asteroids is included. Furthermore, 

• Terrestrial and lunar figure effects, 

• Earth tides, 

• and Lunar physical librations have been considered. 

External non gravitational forces These forces include: 

• Solar radiation and Solar wind pressure. 

• Drag from interplanetary dust. 

Internal (spacecraft) non— gravitational forces These forces include: 

• Thermal radiation. 

• Attitude-control propulsive maneuvers and propellant (gas) leakage from the spacecraft's 
propulsion system. 

• Torques produced from the above forces. 

Model of observation stations An orbit determination has to include a model of the ground 
stations. This is based on: 

• Precession, nutation, sidereal rotation, polar motion, tidal effects, and tectonic plates 
drift. All the informations on tidal deceleration, non-uniformity of rotation. Love num- 
bers, and Chandler wobble have been obtained from LLR, SLR, and VLBI (from ICRF) 
measurements. 

• Model of DSN antennae and their influence on the tracking data. 

Modeling of signal propagation The propagation of the radio signals includes 

• a relativistic model for light propagation including order v^ jc? . 

• and dispersion due to Solar wind and interplanetary dust. 

Codes Four independent codes have been used for the orbit determination: 

• JPL Orbit Determination Program (various generations from 1970 - 2001). 

• The Aerospace Corporation code POEAS (during period 1995 ~ 2001). 

• Goddard Space Flight Center conducted a study in 2003 (data from NSSDC). 

• Code of University of Oslo. 

The definition of these models have to be complemented by a discussion of possible errors. This 
is tantamount to the search for possible conventional explanations of the effect. We present a few 
points only. 

4.3 Discussion of some conventional effects 

In the following we discuss recent and ongoing work on conventional effects which may contribute 
to errors or perhaps also may be responsible for the observed acceleration. Not included here is the 
spin-rotation coupling [22] 



Dust The interplanetary medium consists of (i) interplanetary dust and of (ii) interstellar dust. 
The first one consists of hot wind plasma (mainly protons and electrons) distributed within the 
Kuiper belt (from 30 AU to 100 AU). The density of this plasma has been modeled to be of the 
order pipD < 10^^^ g/cm'^ (Man and Kimura 2000). The interstellar dust which can be distinguished 
from the interplanetary dust by its greater impact velocity has been measured by Ulysses to have a 
density of pisD < 3 • 10"^^ g/cm^. 

The drag acceleration of a spacecraft moving through dust of density p is given by 

fldrag = -KsPVs—^ , (6) 

rris 

where Kg is the satellite's drag coefficient which can be taken to be ~ 2. If we assume the drag 
acceleration to be the observed anomalous acceleration of the Pioneer spacecraft, then this needs a 
density which is 3 • 10^ larger than the interplanetary dust. Therefore dust cannot be the origin of 
the Pioneer acceleration [23]. 

Additional masses in the Solar system Additional masses may not only be present in the 
form of dust but also in form of larger particles. Irrespective of being dust and of the size of these 
particles, any additional mass will act as an additional gravitational field which may decelerate the 
spacecraft when leaving the Solar system. Nieto [24] has calculated analytically the gravitational 
effect of various configurations, that is, shells, thin rings and wedges of various density profiles. He 
obtained that for rings with a density falling off with 1/r as well as a wedge with a density falling 
off like 1/r^ yields a nearly constant acceleration (neglecting discontinuities at the sharp boundaries 
of the matter distributions which are, of course, just results of the mathematical model). However, 
in order the constant acceleration to be of the order of the observed Pioneer acceleration, the mass 
of the thin ring or the wedge has to be of about 100 time the mass of the Earth which is, by far, 
not compatible with the observations of, e.g., comets. 

Accelerated Sun A nongravitational acceleration of the Sun orthogonal to the ecliptic will also 
cause an acceleration toward the Sun. Such an accelerated Sun is the consequence of an exact 
solution of the Einstein equation, the so called C-metric [25]. In the frame of an accelerated Sun, 
the equation of motion for test masses reads 

r + GMQ^ + aQ = 0, (7) 

where r is the distance between the Sun and the test mass. This gives a constant acceleration 
toward the Sun [25]. However, in order obtain an acceleration of the order of the Pioneer anomaly, 
the acceleration of the Sun orthogonal to the ecliptic has to be larger than what would be obtained 
if all radiation of the Sun is emitted in one direction. 

Cosmic expansion Due to the quite good equality apioneer ~ cH where H is the Hubble constant, 
it has been speculated whether the cosmic expansion has some influence on the (i) signal propagation, 
(ii) trajectory of the spacecraft, (iii) the magnitude of the gravitational field inside the Solar system, 
or on (iv) the definition of the distance, that is, the definition of the Astronomical Unit. 

The infiuence of the expansion of the universe on the procedure of Doppler tracking is negligible. 
For an expansion described by an Einstein-de Sitter universe 

ds'^ = -dt^ + R^{t) {dx'^ + dy'^ + dz'^) (8) 



we obtain the conserved quantities Vu{t)R{'t) = const. , where v = v^ = k{u) is the measured 
frequency of a hght ray g{k, k) = 0; u with g{u, u) = 1 is the 4-velocity of an observer at rest in the 
cosmic substrate. To first order in the expansion, 

^u{t) = ^^Mto) ~ (1 - H{t - to)) uuito) , (9) 

where H = R/R is the Hubble constant. This describes the Hubble red shift. 
For massive particles we have the conserved quantity 

R^it)-^ = const. 4» R(t) , ^Vjt) = const. , (10) 

as y/l - V^{t) 



where V, defined by g{u,v) = l/y^l — V'^/c^, is the measured velocity of an object moving with 
4-velocity v along a geodesic, D^v = 0. For small velocities R{t)V{t) = const, implying a slowing 
down of the velocity 

yit2) = ^/ih) = (1 - H{t2 - ti)) V{t,) (11) 

A distance D measured by time-of-flight of light rays is defined hy D = R{r2 — ri) . The measured 
velocity of a moving object is then 

^D = R{r2 - n) + Rr2 = HD + V2 (12) 

where V2 is the velocity of the object measured with respect to the cosmological substrate. As 
a consequence, the trajectory of an object which has constant distance to an observer, = Z) = 
HD + V2 , has to move with velocity 

V2 = -HD (13) 

with respect to the substrate. 

Therefore, the two way Doppler ranging is influenced by three effects: (i) the cosmological 
redshift (9), (ii) the slowing down of the velocity of the spacecraft (11), and the velocity of the unit 
distances with respect to the cosmological substrate (13). Taking all together results in the final 
two-way Doppler ranging signal 

— = -2V{ti){l-H{t2-t{)), (14) 

where ti and ^2 are the cosmological time parameters for the emission and reception of the signal. 
The expansion induced effect is a chirp of Doppler signal related to an acceleration 

a = HV = —cH (15) 

c 

which is by a factor V/c smaller than observed Pioneer acceleration apionccr ~ cH. Therefore, 
Doppler tracking in an expanding universe cannot account for the observed Pioneer Anomaly. 

Furthermore, the influence of the expansion of the universe on the gravitational field of the Sun 
or the planetary orbits is much too small to be of any influence, see Sec. 9. 2. 



Table 1: Data sets analyzed previously and to be analyzed in the future. 



previously used data set 
time span AU distance 



data set to be analyzed 
time span AU distance 



Pioneer 10 
Pioneer 11 



3.1.1987 - 22.7.1998 40 - 70.5 
5.1.1987-1.10.1990 22.4-31.7 



8.9.1973 - 27.4.2002 4.6 - 80.2 
10.4.1974 - 11.10.1994 1.0 - 41.7 



Drift of clocks on Earth Though a drift of clocks by itself is non-conventional physics (the drift 
of atomic clocks with respect to, in our case, gravitational time scales needs some "new physics") it 
may yield a conventional explanation of the Pioneer anomaly. A quadratic drift of the timekeeping 
of clocks on Earth may simulate the Pioneer anomaly [21, 26] 



t 



1 2 

t + 7r«p* 
2c 



(16) 



The numerical value is 2^ap ~ 10~^^ s~^. If we assume such a kind of clock drift then the question 
arises whether this clock drift is consistent with the measurements from other satellites and, in 
particular, with the ranging of satellites. Another question is whether this time drift is also consistent 
with the observations of pulsars and binary systems which also define clocks. Pulsars are very stable 
clocks [27]. Owing to the radiation of gravitational waves the revolution time of binary systems goes 
down. However, also in this case the stability of this process can be defined [28]. A comparison with 
a drift of Earth clocks has not been carried through. 

To sum up, the expansion of our universe seems to be of no relevance for the occurrence of the 
Pioneer acceleration. 



4.4 Outlook 

New activities are planned for the very near future. During the last months the complete set 
of Pioneer data have been recovered [29, 30] and brought into a digital form readable by modern 
computers. These data are now ready for a new data analysis covering all data and the total mission 
duration, see table 1. This new data analysis will be carried through at ZARM and at JPL. It is 
important to find out, e.g., whether the anomalous acceleration was really not present before the 
last flyby. Furthermore, a new Deep Space Gravity Explorer mission has been proposed [31]. 



5 The flyby anomaly 

5.1 The observations 

It has been observed at various occasions that satellites after an Earth swing-by possess a significant 
unexplained velocity increase by a few mm/s. This unexpected and unexplained velocity increase is 
called the flyby anomaly. According to information from [32, 33, 34] the observed flybys are listed 
in Table 2. For the actual data for the Galileo and NEAR flyby see Fig.l. 

The data can be put into diagrams where the velocity increase can be plotted as a function 
of the two orbit parameters eccentricity e and pericentre rp, see Fig. 2. In general, this is a three- 
dimensional plot Av = f{e, Vp). For a plot of this surface four data points are far too less. Therefore, 
in Fig. 2 a plot of the velocity increase Av as function of e and of Vp, respectively, is given. Though 
from four data points it is much too early to draw any serious conclusion one may speculate the 



Table 2: Observed flybys. 



Mission 


agency 


year 


pericentre 


eccentricity 


velocity increase 


Galileo 


NASA 


Dec 1990 


959.9 km 


2.47 


3.92 ± 0.08 mm/sec 


Galileo 


NASA 


Dec 1992 


303.1 km 


2.32 


no reliable data*^ 


NEAR 


NASA 


Jan 1998 


538.8 km 


1.81 


13.46 ± 0.13 mm/sec 


Cassini 


NASA 


Aug 1999 


1173 km 


5.8 


0.11 mm/sec 


Stardust 


NASA 


Jan 2001 


5950 km 




no reliable data*^ 


Rosetta 


ESA 


Mar 2005 


1954 km 


1.327 


1.82 ± 0.05 mm/s 


Hayabusa 


Japan 


May 2004 


3725 km 


?? 


no data available 


MESSENGER 


private 


Aug 2005 


?? 


?? 


no data available'^ 



^ too low orbit with too large atmospheric drag 

thruster activities 
^ US spacecraft operated by a private company 
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(a) Two-way S-band Doppler residuals and range residuals during the first Galileo fiyby 
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(b) Two-way X-band Doppler residuals and range residuals during NEAR flyby 



Figure 1: Galileo and NEAR flyby data (from [33]). 



following: If the velocity increase really is due to an unknown gravitational interaction, then it 
makes sense that (i) the effect goes down with increasing eccentricity, since for larger eccentricity 
the strength and duration of the interaction with the gravitational field of the Earth goes down, 
and (ii) that it also should go down for an eccentricity approaching e = 1 because the transition to 
bound orbits, where no effect has been seen, certainly should show no discontinuity due to many 
kinds of disturbances like drag and non-ideal circumstances like gravitational multipoles, etc. But 
relying on this poor data base this interpretation is pure speculation only. 

The main problem is not just the limited number of flybys for which sufficiently precise data are 
publicly available so that the anomaly can be seen at all. Even these available data suffer from low 
cadence (the anomaly often appears between two data points) and so far only allow an anomaly in 
the speed, but not in the direction of motion etc. to be identified. Precise data at a much higher 
cadence of all the motion parameters of the spacecraft prior to, during and after the flyby would 
allow a qualitatively improved analysis. 



NEAR 




Rosetta 



Figure 2: The velocity increase At; as function of the eccentricity and of the perigee. 



5.2 Error analysis 

As a first remark we give the order of the acceleration leading to the velocity increase. This anoma- 
lous acceleration, estimated by the velocity increase during the time of flight near the Earth, is of 
the order 10~^ m/s^. This is considerably larger than the above discussed Pioneer anomaly. Below 
we will use this acceleration as an approximate value to be compared with disturbing influences. 

It should be kept in mind that the velocity increase has been observed in the two-way Doppler 
measurements as well as in the ranging measurements. 

Before one starts with possible fundamental explanations of this effect, a serious and reliable 
error analysis has to be carried through. This analysis has to cover (i) atmospheric mismodeling, 
(ii) ocean tides, (iii) if the spacecraft becomes charged, then it may experience an additional force 
due to the Earth's magnetic field, (iv) also the interaction of a hypothetical magnetic moment of 
the spacecraft with the Earth's magnetic field may give an additional force, (v) ion plasma drag, 
(vi) Earth albedo, and (vii) Solar wind. Here we give very short first order estimates on these 
various effects which certainly has to be improved, and show that even with very rough pessimistic 
assumptions none of these can be held responsible for the flyby anomaly. 



Atmosphere If a spacecraft of mass m^ and effective area A^ moves with velocity Vg through a 
medium of density p, then it experiences a drag acceleration again given by (6). For a mass 
of 1 t, an area of 2 m^, a velocity of 30 km/s and an atmospheric density at 1000 km height 
of approx p ~ 10^^^ kg/m we get an acceleration of a^rag ^ 4 • 10^^ m/s what is far too 
small to be of any relevance for our problem. Furthermore, this acceleration due to drag has 
the wrong sign. 

Ocean tides The ocean tides will lead to a change of the Earth's surface of the order of 6r it 10 m. 
This means that the corresponding quadrupole part of the Earth's gravitational potential is 
of the order e = 25r/R^ smaller than the monopole part of the Earth, where R(^ is the radius 
of the Earth. Since e ~ 10^^, the corresponding additional acceleration also is factor 10^^ 
smaller than the ordinary acceleration from the monopole part of the Earth's gravitational 
field. The latter being less than 10 m/s , the acceleration due to tides is at most 10^^ m/s 
and, thus, cannot be responsible for the flyby anomaly. 

Solid Earth tides Since Earth solid tides are much smaller than ocean tides, the analysis above 
shows that this cannot cause the effect. 

Charging of the spacecraft In a recent study of charging of the LISA test masses [35] the charg- 
ing has been estimated by 10""'^'^ C. So, for the whole satellite it might be a conservative 
assumption that the charge is less than Q < 10^^ C. A satellite of 1 t carrying a charge Q 
and moving with t; = 30 km/s in the magnetic field of the Earth which is of the order 0.2 G 
will experience an acceleration 10~* m/s far below the observed effect. 

Magnetic moment The force on such a body carrying a magnetic moment m moving in a magnetic 
field B is F = V(m, ■ B). Since the magnetic moment of a spacecraft is not more than 2 A m^ 
and the steepness of the magnetic field can be estimated by \AB/Ax\ < 2 • 10^^ G/m, see 
Fig. 3, the maximum force of a spacecraft is F < A ■ 10^^^ N implying typically a maximum 
acceleration of 4 • 10~^^ m/s which safely can be neglected. 

Earth albedo The Earth albedo causes a pressure on the spacecraft of approx 1 pN/m^ which 
leads, for an effective area of 2 m^ to a force of 2.4 pN. For a mass of the spacecraft of 1 t 
this will give an acceleration of aaibcdo ~ 2.4 • 10~^ m/s what can be neglected compared to 
the searched for effect of 10~^ m/s . 

Solar ■wind The solar wind exerts on spacecraft a pressure of approx 4 /xN/m^ which gives an 
acceleration of max asoiarwind ~ 2.4 • 10~^ m/s which again can be safely neglected. 

Spin— rotation couphng A coupling of the helicity of the radio waves with the rotation of the 
spacecraft and the rotation of the Earth also leads to an effect which simulates a changing 
velocity [22]. This, however, applies to the two-way Doppler data only. Since simultaneously 
also ranging, what is independent of the helicity-rotation coupling, indicated an increase of 
the velocity, spin-rotation cannot be responsible for this observation. 

Also estimates of the influence of the Moon including Moon oblateness, the Sun, other planets, 
relativistic effects, and indirect oblateness of the Earth have been shown to be order of magnitude 
smaller than the observed effect [33]. 

None of these disturbing effects could explain the flyby-anomaly. 
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IGRF-8: result lor 1 orbit: 
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Figure 3: The magnetic field of the Earth as function of attitude. 

5.3 Explanations from "new physics" 

As reported in the paper [33] several non-standard physical models have been used to explain this 
velocity increase. Since these considerations have not been published we just mention the discussed 
models: 

• Non-conservative potential energy. 

• Non-Newtonian gravity (e.g. Moffat, Yukawa, etc). 

• PPN. 

• Modifications of relativity. 

• Torsion, i.e. the e'ps2 model. This model is said to fit the data, but is not compatible with 
stability of planetary orbits. 

As told in [33] none of these could explain the flyby-anomaly. Being so large and so near to the 
Earth, the expansion of the universe should play no role for this anomaly. In Sec. 9 we sketch a 
general approach to a description of the motion of test bodies. Among the terms found there are 
several which may be considered as phenomenological description of a velocity increase. In order to 
be able to pin down a specific interaction term it is necessary to make detailed studies of the nature 
of the velocity increase: 

5.4 Future flybys 

In the near future there will be two fiybys, both by Rosetta [34] 

• Rosetta: flyby on 13 November 2007 (pericentre altitude 4942 km). 

• Rosetta: flyby on 13 November 2009 (pericentre altitude 2483 km). 

We strongly suggest that due to the lack of explanation of the flyby anomaly one should use these 
opportunities in order to carry through a better observation of the Rosetta flybys. A better data 
basis then will enable one to establish a correlation between the observed velocity increase and 
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the orbital parameters like eccentricity, perihelion distance to the Earth, perihelion velocity, or 
inclination. In particular, a continuous observation (Doppler tracking, ranging, positioning, and 
perhaps other data from the spacecraft like temperature, pressure, etc) also should give hints to 
the particular direction of the local acceleration and also on the strength and, thus, to the position 
dependence of the anomalous force. Furthermore, it would also be of great importance to know 
whether after a flyby the direction of the motion of the spacecraft also is different compared with 
the standard result. These informations are extremely important for the search for a conventional 
explanation of this effect or for the modeling of a new force. 

As a consequence, a complete, accurate and continuous observation of spacecraft during flybys is 
very important in order to study the nature of this unexplained velocity increase. It is an advantage 
of these planned observations that they will in no way whatsoever modify the mission scenario. 

Furthermore, since new missions like the Mars Reconnaissance Orbiter will measure very precisely 
the gravitational field of Mars, also the 

• Mars flyby of Rosetta on 25 February 2007 

should be observed as precise as possible in order to see whether also at Mars a velocity increase 
will occur. A Mars flyby would provide an excellent augmentation of the Earth flybys. Since 
Mars possesses other conditions than the Earth (weaker atmosphere, almost no magnetic field, 
other gravitational field, lower thermal radiation, etc.) many competing effects and also possible 
explanation schemes based on the fact that the effect as well as the observation occurs at and on the 
Earth, can be ruled out. Therefore the effect, if it will be observed also at Mars, then will turn out to 
be universal and beyond any doubt and will become an extremely important science case. In the case 
of Mars it is of course impossible to have a continuous observation (at least due to a Mars eclipse) , 
but the initial and final velocity should be measured with the best possible precision. Furthermore, if 
this effect really seems to be existent, then a dedicated mission with a good drag-free control system 
(which can control the acceleration better than 10 m/s for 1 s measurement time) for well-defined 
fiybys at Earth for small perigees and, thus, for extreme parameter values, might be very helpful 
for a even better exploration of this effect. 

6 The increase of the Astronomical Unit 

6.1 The observation 

From the analysis of radiometric measurements of distances between the Earth and the major planets 
including observations from Martian orbiters and landers from 1961 to 2003 a secular increase of 
the Astronomical Unit of approximately 10 m/cy has been reported [36] (see also the article [37] 
and the discussion therein). 

6.2 Search for explanation 

Time dependent gravitational constant and velocity of light This increase cannot be ex- 
plained by a time-dependent gravitational constant G because the G/G needed is larger than the 
restrictions obtained from LLR. 

It has also been speculated that a time-dependent change in the velocity of light can be respon- 
sible for this effect. Indeed, if the speed of light becomes smaller, than ranging will simulate a drift 
of distances. However, a inspection of Kepler's third law 

^ = ^ (17) 
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shows that, if one replaces the distance a by a ranging time a = ct, then effectively the quotient G/c^ 
appears. Only this combination of the gravitational constant and the speed of light governs the ratio 
between the orbit time, in our case the orbit time of the Earth. Consequently, a time-dependent 
speed of light is equivalent to a time-dependent gravitational constant. Since the latter has been 
ruled out to be possibly responsible for an increase of the Astronomical Unit, also a time-dependent 
speed of light has to be ruled out. 

Cosmic expansion The influence of cosmic expansion by many orders of magnitude too small, see 
Sec. 9. 2. Neither the modification of the gravitational field of the Sun nor the drag of the planetary 
orbits due to the expansion is big enough to explain this drift. 

Clock drift An increase of ranged distances might also be due to a drift of the time scale of the 
form t ^ t + at^ for a > 0. This is of the same form as the time drift needed to account for the 
Pioneer anomaly. From Kepler's third law one may ask which a is suitable in order to simulate the 
increase of the Astronomical Unit. One obtains a ~ 3 • lO"^'^ s~^ what is astonishing close to the 
clock drift needed for a clock drift simulation of the pioneer anomaly, see Eq.(16) and below. 

7 The quadrupole and octupule anomaly 

Recently an anomalous behavior of the low-/ contributions to the cosmic microwave background has 
been reported. It has been shown that (i) there exists an alignment between the quadrupole and 
octupole with > 99.87% C.L. [38], and (ii) that the quadrupole and octupole are aligned to Solar 
system ecliptic to > 99% C.L. [39]. No correlation with the galactic plane has been found. 

The reason for this is totally unclear. One may speculate that an unknown gravitational field 
within the Solar system slightly redirects the incoming cosmic microwave radiation (in the similar 
way as a motion with a certain velocity with respect to the rest frame of the cosmological background 
redirects the cosmic background radiation and leads to modifications of the dipole and quadrupole 
parts). Such a redirection should be more pronounced for low-/ components of the radiation. It 
should be possible to calculate the gravitational field needed for such a redirection and then to 
compare that with the observational data of the Solar system and the other observed anomalies. 

8 Summary of anomalies 

8.1 Summary 

In our opinion, these above described anomalies split into three groups according to their observa- 
tional status: (1) The observations related to dark matter and dark energy are beyond any doubt, 
(2) the Pioneer anomaly and the flyby anomaly are on a good basis while (3) the increase of he 
Astronomical Unit as well as the quadrupole and octupole anomaly are still under debate. 

This list of anomalies and unexplained phenomena immediately induces a bunch of tasks and 
questions: 

• For each of these phenomena, except for the Dark Energy and Dark Matter, one still should 
try to find a systematic cause. 

• One also should try to find conventional explanations of these effects within General Relativity 
leading perhaps to an effect within conventional physics not considered up to now. 

• A big step in understanding these effects might be to find a relation between two or more of 
the unexplained phenomena (all of these anomalies very probably are no isolated phenomena) . 
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Table 3: List of anomalies and their status 



Anomaly 


Observational status 


Interpretation 


Dark Energy, 
Dark Matter 


well established 


under discussion 


Pioneer Anomaly, 
Flyby anomaly 


quite well established 


unclear 


Quadrupole anomaly, 
increase of AU 


unclear 


unclear 



• This is tantamount to the important question whether all or at least some of these effects have 
a common cause. 

• Motivated by the surprising almost perfect equality of apioneer = cH one should analyze the 
influence of the cosmic expansion on the physics within gravitationally bound systems. There 
already appeared quite some literature related this topic [40, 41, 42, 43], and references therein. 

• Since there seem to exist some strange phenomena within our Solar system, one should look 
whether it might be possible to observe similar things in other gravitating systems like binary 
systems, binary pulsars, stars moving around the black hole in the center of our galaxy, etc. 

• Is something wrong with weak gravity or gravity at large distances? 

• Furthermore, one should propose new dedicated missions and space experiments. 

8.2 Other anomalies? 

There is one further observation which status is rather unclear bit which perhaps may fit into the 
other observations. This is the observation of the return time of comets: Comets usually come 
back a few days before they are expected when applying ordinary equations of motion. The delay 
usually is assigned to the out gassing of these objects. In fact, the delay is used for an estimate of the 
strength of this outgassing. On the other hand, it has been calculated in [44] that the assumption 
that starting with 20 AU there is an additional acceleration of the order of the Pioneer anomaly 
also leads to the effect that comets come back a few days earlier. It is not clear whether this is a 
serious indications but a further study of the trajectories of comets certainly is worthwhile. 

9 Ways to describe the effects 

Many approaches have been attempted to explain these anomalies. In most of these attempts a 
link between one phenomenon and the issues of (i) the influence of the expansion of the universe 
on the physics within our Solar system, (ii) dark energy, and (iii) dark matter has been tried to 
establish. We like to emphasize again that indeed it should be the first thing to explore whether 
there are links between these various observed and unexplained phenomena. It should be strange if 
all these unexplained phenomena will be really independent of each other, that is, are not linked by 
a common (perhaps new) physical principle. 

However, it also seems to be a rather difficult task to find similarities between the Pioneer 
anomaly and the flyby anomaly. 
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9.1 Cosmological constant 

One first very simple attempt to generalize General Relativity and to incorporate also findings 
related to dark energy is to describe the physics of gravitating bodies within a theory involving a 
cosmological constant. The corresponding fields equations are 

Rtiu - ^gfiuR + g^ivJ^ = nT^u (18) 

where R^^ is the Ricci tensor, R the Ricci scalar, and T^^y the energy momentum tensor. 

In a first step one may use this equation in order to describe the gravitational field of the Sun 
and to describe the physics within the Solar system (Perihelion shift, light deflection, gravitational 
redshift, gravitational time delay, geodetic precession) in [45] as well as by lorio [46] and Sereno 
and Jetzer [47]. It comes out that the ordinary cosmological constant A = Aq = 10~^^ m~^ has no 
observable influence on all the Solar system tests. One also can show that if one assumes a constant 
A such large that it may explain the Pioneer acceleration, then this is ruled out by the perihelion 
shift observation [45]. Therefore, the cosmological constant cannot be an explanation for the Pioneer 
anomaly. In the paper by Jetzer and Sereno [48] also the influence of the cosmological constant on 
the motion of binary systems has been evaluated with the result that in the future binary systems 
may be precise enough to "see" the cosmological constant. Some of these effects have also been 
considered in a Kerr-de Sitter spacetime [49]. 

The calculation of Solar system effects in a Schwarschild-de Sitter space-time is certainly only a 
simple first step. One may redo this kind of calculation in a more general context like quintessence 
[50, 15] (for a first step in this direction see [51]), in varying G scenarios [52], in dilaton scenarios 
[53, 54], and in braneworld models [55]. 

9.2 The influence of the expansion of the universe 

We show that the expansion of the universe hat no measurable influence whatsoever on the physics 
within the Solar system. This includes the modifications of the gravitational field created by the 
Sun and the planetary orbits. 

Modified gravitational field of the Sun It is also conceivable that the cosmic expansion may 
weaken the gravitational field of the Sun. However, it can be shown [42] that the corresponding 
effect is far beyond being observable. Starting from an expansion of the local metric around the 
cosmological background metric b^i, 

OfMu = bf^u + hf^u , where hf^^ <■ b^u ■ (19) 

we obtain linearized Einstein equations for h^i, [56] 

gf^DpD^h^, + 2(7^'^i?%p,/i,, = WttGT^, , (20) 

The static solution for a small spherically symmetric mass distribution is given by 

2GMcos(\/6|^|r) GM , 9. .9 ^ 

At lowest order we obtain the standard Newtonian potential with the measured distance R{t)r. 
Since the modification is quadratic in the Hubble constant, the Newtonian potential practically 
does not participate in the cosmic expansion. This confirms the findings in [36]. 
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Modified planetary orbits If the planetary orbits may expand due to a drag from the cosmic 
expansion, then this might be interpreted as if the length unit, Astronomical Unit, increases so that 
it might appear that other distances may become smaller. However, from the theorem of adiabatic 
invariance [57] it can be shown that the parameters of the planetary orbits remain the same with 
very high stability [42]. This stability is characterized by the factor exp(— t/T) where t is the time 
scale of the planetary orbit and T the time scale of the temporal change of the environment, in our 
case the Hubble time. Therefore, t/T is of the order lO^'^ so that the planetary orbits are extremely 
stable. 

Another approach is based on the geodesic deviation which relates the relative acceleration of 
freely moving particles with the curvature of space-time, in our case with the expansion of the 
universe. This has been taken in [58] as the basis for estimates that this deviation is much too small 
to be of any influence in the Solar System or even within galaxies. For a broader discussion of this 
approach, see the recent valuable review of Carrera and Giulini [43]. 

9.3 General approach to describe a modified particle dynamics 

Gravity describes or, equivalently, can be characterized by the behavior of light and point particles 
[59]. Light is related to a space-time metric while point particles are related to a geodesic equation. 
In Riemannian geometry, the mathematical model for General Relativity, the equation of motion 
for particles and the behavior of light are related: The equation of motion for particles is completely 
determined by the space-time metric. For general analysis of anomalous gravitational effects it 
therefore may be appropriate to start a by a general ansatz for the metric and of the equation of 
motion for particles^. 

In a weak field approximation the metric can be written in the form 

where U, V, and goi = gio are assumed to be small quantities. U may be identified with the 
Newtonian potential. We denote ds'^ = g^^dx'^dx'^ . 

By definition, the readout of moving clocks is given by 

T= f ds ^ f{l-U + x'^ + Vx^)dt+ f h-dx (23) 

^worldlinc of clock J J 

which is the proper time of the clock. A clock at rest will measure the gravitational redshift given 

by 

- = \lf^,-' + U{x,)-U[x,). (24) 

Now we are going to give a general description of the equation of motion for particles. This 
approach is according to the philosophy of axiomatic approaches to theories of gravity, see, e.g., 
[59, 60]. The equation of motion for a particle respecting the Universality of Free Fall can be 
represented as 

= v''d^v'' = H^{x,v) = {p^a}v^'v^ + -/^'{x,v) (25) 

^In the axiomatic approach [59] the mathematical relation between the point particle dynamics and the properties of 
light came in through a compatibility condition. One feature of this condition is a causal dynamics of point particles; 
particles are not allowed to move faster than light. The requirement of this condition implies a rather restricted 
geometrical structure, namely a Weylian structure. Since in [59] this condition has been applied to autoparallel curves 
and since we are here more general in admitting also non-linear connections as they appear, e.g., in a Finslerian 
context, the resulting allowed particle dynamics are certainly still more general than a Weylian structure. However, 
this has to be worked out explicitly. 
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where {p^} is the ChristofFel connection and 7^(x,f) some vector valued function of the position 
and the velocity. Since no particle parameter enters this equation, the Universality of Free Fall 
automatically is preserved. 

From this general equation of motion, we can derive the 3~acceleration 

d?x^ / r ,■ -, ( n -^ dx^ \ dx^ dx^ 1 / ,• , , dx* 



dt"^ 



/(. . . . . dx''\ dxi^ dx" 1 f i dx' Q 

\ds ) 

{«U-{o»o}^)-2({,,U-{,«^)^-({A}-{,M^) 



dt j dt dt 

« diU + {dihj - djhi)x^ +x^diV + xV + T^ + T\x^ + Tkx^x^ + . . . , (26) 

^~^ ^ ^ ' 

Newton Lense-Thirring 

where we neglected all relativistic corrections since these play no role in the Pioneer and flyby 
anomalies. Note also that the Universality of Free Fall is respected. However, it is no longer 
possible to make a transformation to a coordinate system so that gravity disappears at one point 
(Einstin's elevator is not possible). Physically this means that, e.g., the acceleration of a body 
toward the Earth can depend on the velocity of the body: Differently moving bodies feel a different 
gravitational acceleration which, however, does not depend on the composition or the weight of the 
(point-like) body. 

The first term in (26) is the ordinary Newtonian acceleration and the second term the action of 
the gravitomagnetic field on the orbit of a satellite which has been observed by LAGEOS with a 10% 
accuracy [61]. This field also acts on spinning objects like gyroscopes and should be confirmed by 
GP-B with an accuracy better than 1%. The other terms are hypothetical terms beyond ordinary 
post-Newtonian approximation. 

The V term which can be motivated by a running coupling constant to be proportional to the 
distance, y ~ r^ has been introduced by Jaeckel and Reynaud [62] in order to describe the constant 
anomalous Pioneer acceleration. The other terms, most of then are velocity dependent, have not yet 
been analyzed. The influence of an arbitrary force on the trajectories of planets has been analyzed 
recently in [63] with the main conclusion that any radial force which might be considered as being 
responsible for the Pioneer anomaly is not compatible with the recent analysis of the motion of 
the outer planets. This indicates that the modification of the equation of motion should include 
velocity-dependent terms. 

The coefficients T*^ depend on the position only and may vanish for vanishing gravitating 
mass. Therefore, the coefficients can contain M, r, r*, and di only. Accordingly, these coefficients 
can be of the form 

T} = A21—, ^ + A22^6) (28) 

: . GMrW^r^ ^ GM r\ , GM r\. 
T!.fc = A^i— 3— +^32^ Sjk + A^i^— 51. (29) 

1/ fv*^ fy*'-J fy*^ fy* '■' ty* ^ ly* 

Here it is understood that the infiuence of the gravitating body, that is, all the T-coefficients, 
vanishes at spatial infinity. This is certainly true for a description of the flyby anomaly but may be 
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relaxed for the Pioneer anomaly. Terms which do not vanish at spatial infinity but are of Newtonian 
form at small distances are polynomials r', / > 1. 
The above terms lead to accelerations 

. GMr^ , , 

■£ = Au^ (30) 

, GMr^r-r ^ GM r' 
x' = ^21^5 ^ + ^22 



, , CM rW-f ^ GM f\ , , 

= (^21+^22)^:^^^ +^22^:^^ (31) 

, GMrHr-rf , GM r^ o , GM r\ , , , 

where r^^ = r^ — r'^ir ■ r)/r'^ is the component of the body's velocity orthogonal to the connecting 
vector r, and the Aij are some numerical factors. 

The first term associated with An is of Newtonian form and can be combined with the already 
existing one which amounts to a redefinition of the gravitational constant. The A22 term describes 
an additional acceleration in direction of the velocity. It fades away for large r. The A21 term 
projects the component of the velocity which is parallel to the connecting vector and leads to an 
acceleration in direction of the connecting vector. This term vanishes at the perigee. 

These A2i~ and ^22^terms may be chosen in such a way that they have the potential to describe 

an increase of the velocity during a fiyby: Near the Earth the ^22^term is dominant since there 

the connecting vector is more or less orthogonal to the velocity vector. For large r both terms 

contribute. That is 

^ GMe 

A22 -^ for r RS rperigee 

(^21+ ^22) — o for r large 

Therefore, in principle it is possible to have an acceleration near the perigee (for A22 > 0) and a 
deceleration for large distances (for A21 + A22 < 0). Note that for a typical perigee and velocity at 
perigee the acceleration at perigee for A22 = 1 is about 10"'^ m/s what is just the value given for 
a typical Earth flyby, see Sec. 5. 2. However, this model does not include the Pioneer deceleration 
because the acceleration is not constant for large r. What is needed now is a general discussion of 
the influence of a term of the form (31) on general features planetary and satellite orbits, e.g., the 
perihelion shift, and to compare this with observations. Corresponding work is in progress. 

The A32 term just adds to the iP'diV term. The A^i and the ^33 term both first project the 
velocity in direction of the connecting vector and then make out of this an acceleration in direction 
of the connecting vector and in direction of the velocity. These Asj-terms are about 5 orders of 
magnitude smaller than the yl2i-terms and can, thus, play no role in an explanation of the flyby 
anomaly. The higher order terms will be of more complicated but similar structure. 

In general, this equation of motion does not respect energy conservation: multiplication of (26) 
with the velocity yields 



-<\x^-U) = 2x^V + x-r + Tix^i^ + TijjPi^i^ + ... (34) 



dt 



Therefore, the terms on the right hand side might be candidates for effects reducing or enlarging 
the kinetic energy of moving bodies and, thus, may play a role in the description of the fiyby or the 
Pioneer anomaly. 
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It should be clear from the independence of the metric from the equation of motion for point 
particles, that it is necessary both to track position and velocity of the satellite and to have a clock 
on board in order to determine all components of the space-time metric. 

10 Summary and outlook 

We collected the anomalies related to the physics of the Solar system and discussed to some ex- 
tend the error sources and possibilities to explain these anomalies. In particular, we tried to find 
similarities or fundamental differences between these anomalies. 

As final statement we like to stress that there are at least three important science cases related 
to the exploration of these anomalies which we strongly suggest to be tackled in the near future: 

1. Analysis of the complete set of Pioneer data. 

2. Continuous and complete (velocity, distance, time on board, and direction) observations of 
future flybys. 

3. Search for clock drifts by comparison of clock rates on Earth with clocks defined by astrophys- 
ical systems. 
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